Magnetism and superconductivity 
in quark matter 

Toshitaka Tatsumi", Eiji Nakano^ and Kanabu Nawa° 



o 
o 

(N 

P3 ' ^ Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 



February 2, 2008 

> 
(N 

^ , Abstract 

' Magnetic properties of quark matter and its relation to the mi- 

• croscopic origin of the magnetic field observed in compact stars are 

I studied. Spontaneous spin polarization appears in high-density region 

due to the Fock exchange term, which may provide a scenario for the 
behaviors of magnetars. On the other hand, quark matter becomes 
C^. unstable to form spin density wave in the moderate density region, 

^ I where restoration of chiral symmetry plays an important role. Coex- 

istence of magnetism and color superconductivity is also discussed. 

X, 

1 Introduction 

QCD has been believed to be the basic theory of strong interaction and 
there are many successful consequences about the properties of hadrons and 
their interactions. Recently many studies have been devoted to figure out 
the phase diagram of QCD in temperature (T) - density (pb) plane [HE]- 
At high temperature or high density, quarks confined inside hadrons should 
be liberated to form matter consisting of quarks and gluons (deconfinement 
transition). Such basic constituents, especially quarks exhibit interesting 
properties there as electrons in condensed matter through many-body dy- 
namics; one of the interesting possibility is phase transition as temperature 
or density changes. When we emphasize the low T and high ps region, the 



1 



subjects are sometimes called high-density QCD. The main aims in this field 
should be to elucidate the new phases and their properties, and to extract 
their symmetry breaking pattern and low-energy excitation modes there on 
the basis of QCD. On the other hand, these studies have phenomenologi- 
cal implications on relativistic heavy-ion collisions and compact stars like 
neutron stars or quark stars 

Color superconductivity (CSC) should be very popular [Hll]. Its mecha- 
nism is similar to the BCS theory for the electron-phonon system in which 
the attractive interaction of electrons is provided by phonon exchange and 
causes the Cooper instability near the Fermi surface. As for quark matter, the 
quark-quark interaction is mediated by colored gluons, and is often approxi- 
mated by some effective interactions, e.g., the one-gluon-exchange (OGE) or 
the instanton-induced interaction, both of which give rise to the attractive 
quark-quark interaction in the color anti-symmetric 3 channel. Many people 
believe that it is robust due to the Cooper instabihty even for small attractive 
quark-quark interaction in color 3 channel. 

Here we'd like to address another interesting property of quark matter, 
magnetic properties of quark matter. We shall see various types of magnetic 
ordering may be expected in quark matter at finite density or temperature. 
They arise due to the quark particle-hole {p — h) correlations in the pseudo- 
scalar or axial- vector channel. 

Phenomenologically the concept of magnetism should be directly related 
to the origin of strong magnetic field observed in compact stars jB]; e.g., it 
amounts to O(IO^^G) at the surface of radio pulsars. Recently a new class 
of pulsars called magnetars has been discovered with super strong magnetic 
field, Bs ~ lO^^^^^G, estimated from the P — P curve ff^. First observations 
are indirect evidences for the superstrong magnetic field, but discoveries of 
some absorption lines stemming from the cyclotron frequency of protons 
have been currently reported |S|; when it is confirmed that they originate 
from protons, they give a direct evidence for the superstrong magnetic field. 

Table 1: Surface magnetic field and the radius of stars by the conservation 
of the magnetic flux. 

~ Bs[G\ R[cm \ 

Sun (obs.) 10^ 10^° 
Neutron star 10^^ 10^ 
Magnetar 10^^ 10^ 
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The origin of the strong magnetic field in compact stars has been a long 
standing problem since the first discovery of a pulsar A simple working 
hypothesis is the conservation of the magnetic fiux and its squeezing during 
the evolution from a main-sequence progenitor star to a compact star; B oc 
with R being the radius. Taking the sun as a typical main-sequence 
star, we have B ~ lO^G and R ~ lO^^cm. If it is squeezed to a typical 
radius of usual neutron stars, R ~ 10km, the conservation of the magnetic 
fiux gives 10^^ G, which is consistent with the observations for radio pulsars. 
However, we find R ~ 100m to explain B ~ lO^^G observed for magnetars, 
which may lead to a contradiction since the Schwatzschild radius is O(lkm) 
for the canonical mass of 0{Mq), which is much larger than R. 

Since dense hadronic matter should be widely developed inside compact 
stars, it would be reasonable to inquire a microscopic origin of such strong 
magnetic field: ferromagnetism (FM) or spin polarization is one of the can- 
didates to explain it. Makishima also suggested the hadronic origin of the 
magnetic field observed in binary X-ray pulsars or radio pulsarsjH], since it 
looks no field decay in these objects. 

When we consider the magnetic-interaction energy by a simple formula, 
-Emag = fJ'iB with the magnetic moment, fii = ej/(2mj), we can easily es- 
timate it for B = O(IO^^G) (see Table 2); it amounts to several MeV for 
electrons, while several keV for nucleons and 10 keV- IMeV for quarks de- 
pendent on their mass. This simple consideration may imply that strong 

Table 2: Magnetic interaction energies £'mag for lO^^G and the typical energy 
scales Etyp in electron, nucleon and quark systems. 





electron 


proton 


quark 


mi [MeV] 


0.5 


10^ 


1- 100 


^mag[MeV] 


5 - 6 


2.5 X 10-3 


2.5 X 10-2 - ; 


-Etyp 


keV 


MeV 


MeV 



interaction gives a feasible origin for the strong magnetic field, since its typi- 
cal energy scale is MeV. The possibility of ferromagnetism in nuclear matter 
has been elaborately studied since the first discovery of pulsars, but negative 
results have been reported so far 10 . We consider here its possibility in 
quark matter as an alternative in light of recent development of high-density 
QCD lllj. 

If FM is realized in quark matter, there should be some interplay with 
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CSC; we examine a possibility of the coexistence of FM and CSC in quark 
matter, where we shall see an interplay between particle-particle and particle- 
hole correlations. As far as we know, interplay between the color supercon- 
ducting phase and other phases characterized by the non-vanishing mean 
fields of the spinor bilinears has not been explored except for the case of 
chiral symmetry breaking [T^ . 

It would be worth mentioning in this context that ferromagnetism (or spin 
polarization) and superconductivity are fundamental concepts in condensed 
matter physics, and their coexistent phase has been discussed and expected 
for a long time [12]. As a recent progress, superconducting phases have been 
discovered in some ferromagnetic materials and many efforts have been made 
to understand the coexisting mechanism ^Ij. In this phenomenon itinerant 
electrons may be responsible, while its mechanism is not fully elucidated yet. 
Many people believe that the electron Cooper pair should be P— wave, since 
this type can be compatible with spin polarization. We can easily consider 
the similar situation in quark matter. Since the volumes of the Fermi seas of 
quarks with different spins result in being different due to the net presence of 
magnetization, we could not construct a quark Cooper pair in a usual manner 
as = 0^. Instead, we consider the = 0~ pairing with orbital angular 
momentum L = 1 and total spin S = 1. For the 5* = 1 state we further 
consider two possibilities: spin-parallel pair or spin-anti-parallel pair. We 
first discuss the former case in detail, which may have a direct resemblance 
to the electron case. Subsequently we briefly sketch our idea about the 
former case. Anyway we shall see the gap functions become anisotropic in 
the momentum space like in ^He or nuclear matter |^ 1^ . 

We discuss another magnetic aspect in quark matter at moderate den- 
sities, where the QCD interaction is still strong and some non-perturbative 
effects still remain. One of the most important phenomena observed there 
is restoration of chiral symmetry. In the vacuum chiral symmetry is spon- 
taneously broken to give finite mass for quarks or nucleon; we may bear in 
mind such a picture that the vacuum is in a kind of superconducting phase 
with massless quark (g)-anti-quark (g) pair condensate, and the gap opened 
at the top of the Dirac sea corresponds to the finite mass. As a consequence 
the vacuum does not possess chiral symmetry any more. At finite densities 
the suppression of qq excitations due to the existence of the Fermi sea gives 
rise to restoration of chiral symmetry at a certain density, and many people 
believe that deconfinement transition occurs at almost the same density. 

There have been proposed various types of the p-h condensations at mod- 
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erate densities ^Hl QE]; in which the p-h pair in scalar or tensor channel 
has the finite total momentum indicating standing waves (the chiral density 
waves). The instability for the density wave in quark matter was first dis- 
cussed by Deryagin et al. [T3] at asymptotically high densities where the 
interaction is very weak, and they concluded that the density-wave instabil- 
ity prevails over the BCS one in the large N^. (the number of colors) limit 
due to the dynamical suppression of colored BCS pairings. 

In general, density waves are favored in 1-D (one spatial dimension) sys- 
tems and have the wave number Q = 2kp according to the Peierls instability 
[TTl HH], e.g., charge density waves (CDW) in quasi-l-D metals jTH]. The 
essence of its mechanism is the nesting of Fermi surfaces and the level repul- 
sion (crossing) of single particle spectra due to the interaction for the finite 
wave number. Thus the low dimensionality has a essential role to produce 
the density-wave states. In the higher dimensional systems, however, the 
transitions occur provided the interaction of a corresponding {p-h) channel 
is strong enough. For the 3-D electron gas, it was shown by Overhauser 
j2ni 121] that paramagnetic state is unstable with respect to the formation of 
the static spin density wave (SDW), in which spectra of up- and down-spin 
states are deformed to bring about the level crossing due to the Fock ex- 
change interactions, while the wave number does not precisely coincide with 
2kp because of the incomplete nesting in higher dimension. 

We shall see a kind of spin density wave develops there, in analogy with 
SDW mentioned above.. It occurs along with the chiral condensation and is 
represented by a dual standing wave in scalar and pseudo-scalar condensates 
(we have called it 'dual chiral-density wave', DCDW). DCDW has different 
features in comparison with the previously discussed chiral density waves [13 
ITB] . One outstanding feature concerns its magnetic aspect; DCDW induces 
spin density wave. 

2 Ferromagnetism in QCD 
2.1 A heuristic argument 

Quark matter bears some resemblance to electron gas interacting with the 
Coulomb potential; the one gluon exchange (OGE) interaction in QCD has 
some resemblance to the Coulomb interaction in QED , and color neutrality 
of quark matter corresponds to total charge neutrality of electron gas under 
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the background of positively charged ions. 

It was Bloch who first suggested a mechanism leading to ferromagnetism 
of itinerant electrons within the Hartree-Fock approximation |22j . The mech- 
anism looks very simple but largely reflects the Fermion nature of electrons 
in a model-independent way. Since there works no direct interaction be- 
tween charged particles as a whole, the Fock exchange interaction gives a 
leading contribution. Then it is immediately conceivable that a most at- 
tractive channel is the parallel spin pair, whereas the anti-parallel pair gives 
null contribution (see Eq. (jH)) below). This is nothing but a consequence of 
the Pauli exclusion principle: electrons with the same spin cannot closely 
approach to each other, which efficiently avoid the Coulomb repulsion. Thus 
a completely polarized state is favored by the interaction. On the other hand 
a polarized state should have a larger kinetic energy by rearranging the two 
Fermi spheres. Thus there is a trade-off between the kinetic and interaction 
energies, which leads to a spontaneous spin polarization (SSP) or FM at a 
certain density. Subsequently it has been proved that Bloch's idea is qualita- 
tively justified, but the critical density can not be reliably estimated without 
examining the higher-order correlation diagrams [2211211; especially the ring 
diagrams have been known to be important in the calculation of the sus- 
ceptibility of electron gas. Recently the possibility of ferromagnetism in the 
electron gas has been studied by the quantum Monte Carlo simulation and 
it has been shown that the electron gas is in ferromagnetic phase at very low 
electron density |2S1- Authors in ref. [2ni have confirmed it experimentally. 

One of the essential points we learned here is that we need no spin- 
dependent interaction in the original Lagrangian to see SSP: a symmetry 
principle gives rise to a spin dependent interaction. 

Then it might be natural to ask how about in QCD. We list here some 
features of QCD related to this subject. (1) the quark-gluon interaction in 
QCD is rather simple, compared with the nuclear force; it is a gauge interac- 
tion like in QED. (2) quark matter should be a color neutral system and only 
the Fock exchange interaction is also relevant like in the electron system. (3) 
there is an additional flavor degree of freedom in quark matter; gluon ex- 
change never change flavor but it becomes effective through the generalized 
Pauli principle. (4) quarks should be treated relativistically, different from 
the electron system. 

The last feature requires a new definition and formulation of SSP or FM 
in relativistic systems since "spin" is no more a good quantum number for 
relativistic particles; spin couples with momentum and its direction changes 
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during the motion. It is well known that the Pauli-Lubanski vector is 
the four vector to represent the spin degree of freedom in a covariant form; 
the spinor of the free Dirac equation is the eigenstate of the operator, 

W-a = -^75#, (1) 
where a 4-axial- vector is orthogonal to k s.t. 

„« = !^ (2) 



m{Ek + m) ' m 

with the axial vector ^. We can see that a'^ is reduced to a three vector (0, C) 
in the rest frame, where we can allocate <^ = (0, 0, ±1) to spin "up" and 
"down" states. Thus we can still use <^ to specify the two intrinsic polarized 
states even in the general Lorentz frame. To characterize the degeneracy of 
the plane wave solution u^°'\k) (a = 1, 2) for a positive energy state, we can 
use such spinors u^°'\k) that are eigenstates of the operator —W ■ a/niq. for 
the standard representation of u^°'\k) p7j . 

- = ±i«'"'(fc). (3) 

llbq Z 

Accordingly, the polarization density matrix p(/c, C,) is given by the expres- 
sion, 

p(fc,C) = ^(^ + m,)P(a), P(a) = 1(1 + 75^), (4) 

which is normalized by the condition, tTp{k, = 1 |28j . 

Consider the spin-polarized quark liquid with the total number density 
of quarks ^ ; we denote the number densities of quarks with spin up and 
down by n+ and ?t,_, respectively , and introduce the polarization parameter 
p by the equations, n± = \ng{l ±p), under the condition < p < 1. We 
assume as usual that three color states are occupied to be neutral for each 
momentum and spin state. The Fermi momenta in the spin-polarized quark 
matter are then kp = kp{l ±pY^^ with kp = (Tr^rig)-'^/^. The kinetic energy 
density is given by the standard formula. 



167r2 

i=± 



k'pE'p{2k'^ + ml)-mUn' ^ 



(5) 



"'^We, hereafter, consider one flavor quark matter, since the OGE interaction never 
changes flavors. 
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with the Fermi energy Ep = (m^ + kpY^'^. 

Let us consider the OGE interaction between two quarks with momenta, 
k and q, and spin vectors, and <^', respectively. The color symmetric matrix 
element A^^c qC given only by the exchange term; the direct term vanishes 
because the color symmetric combinations (~ trAa) does not couple to gluons. 
Thus 

-^kCqC = -/^tr(A„/2A,/2)n(^')(q)7,n(«(p)n(«(p)7^«(^')(q) 



= l9'\tTMkX)rpiq,C)]j^^„ (6) 

by the use of Eq. Q. If we choose both ^ and C in parallel along the z axis, 
C = C = (0,0, ±1), we have the spin-nonflip amplitude 7V1*'^°"/'*p^ while if 
we choose them in anti-parallel, ^ = — C', we have the spin-flip amplitude 
_/\/js,//ip EaQ]^ form of the spin-nonflip or spin-flip amplitude is complicated, 
but their average gives a simple form. 



which is nothing but the matrix element for the unpolarized case In 
the nonrelativistic limit, ruq ^ |p|, |q|, the matrix element is reduced to the 
form, 

■^kcqc - -^9 — TT, — — ' V^) 



so that there is no correlation between quarks with different spins. On the 
other hand, there is some correlation included in the relativistic case. 

After summing up over the color degree of freedom and performing the 
integrals of the color symmetric matrix element A^^c qC '-'^^^ Fermi seas of 
spin up and down quarks, we have the exchange energy density consisting 
of two contributions, 

_ nonflip I flip (q\ 
'^ex -r Cgj,. . i^^j 

In the nonrelativistic case, the spin-flip contribution becomes tiny and 
the dominant contribution for the OGE energy density in Eq.ijH)) comes from 
the spin-nonflip contribution (see Eq. (jH))), 

^e.~-^{(l+#/^ + (l-|>r/^}. (10) 
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The exchange energy is negative and takes a minimum at p = 1. The form 
of the energy density (llUj) is exactly the same as in electron gas. It is the 
difference of density dependence between the contributions given in Eqs. (0) 
and (jlOp which causes a ferromagnetic instability; this mechanism was first 
pointed out by Bloch for electron gas ^2.. 

In the relativistic case there are some different features from the nonrela- 
tivistic case. First, there is a spin-flip contribution due to the lower compo- 
nent of the Dirac spinor even for the Coulomb-like interaction. Secondly, the 
transverse (magnetic) gluons becomes important, where the spin-flip effect is 
prominent. Finally, the density dependence of kinetic energy as well as the 
exchange energy is very different Before discussing the general case, we 
consider the relativistic limit, kp S> m^; the Fock exchange-energy density 
looks like 

6e. - ^4 {(1 + pY^' + (1 - P)'^' + 2(1 - pY'} , (11) 

which is a decreasing function and takes a minimum again at p = 1. This is 
due to the characteristic feature of the spin-flip and spin-nonflip interactions: 
both give a repulsive contribution in the relativistic limit and there is no 
spin-flip interaction in the polarized state {p = 1). Thus, ferromagnetism 
in the relativistic limit arises by a different mechanism from that in the 
nonrelativistic case. 

In Fig.l a typical shape of the total energy density, etot = ^kin + ^ex, is 
depicted as a function of the polarisation parameter p, e.g. for the parameter 
set ,mq = SOOMeV of the s quark and ac = 2.2 as in the MIT bag model [^2] 
^. We can see that paramagnetic quark matter {p = 0) becomes unastable 
as density decreases, and ferromagnetic phase is favored at a certain density 
between 0.1 and 0.2 fm~^. This phase transition is of weakly first-order and 
the completely polarized {p = 1) state appears at the critical density. 

To figure out the features of the ferromagnetic transition, we study other 
quantities. For small p <^ 1, the energy density behaves like 

etot - etotip = 0) = x''p' + 0{p^) (12) 

with = Xkin + Xex- X is proportional to the magnetic susceptibility, and 
its sign change indicates a ferromagnetic transition, if it is of the second order. 

^The difficulties to determine the values of these parameters have been discussed in ref. 
|33| . and we must allow some range for them. 
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Figure 1: First-order phase transition Figure 2: Critical lines as functions of 
to the ferromagnetic state. quark-number density. 

It consists of two contributions: the kinetic energy gives Xkin ~ ^f/IStt^-E^) 
(c.f. the Pauli paramagnetism), which changes from XkL ~ 0{k%) at low 
densities to Xkln ~ 0{k%) at high densities. On the other hand, the Fock 
exchange energy gives 



^ex - g^3 



^ ^2 _ - ^^l^P 1^ f EF±kF \ ^ ^mgkl 



El E% \ m, ) ?,E%{EF + m, 



(13) 



pUj . which is reduced to 



x-i ~ -^4 (14) 

in the nonrelativistic limit, Pf ^q- In the relativistic limit, Pf ^ f^q, it 
behaves like 

where the first term stems from the spin-nonflip contribution, while the sec- 
ond term from the spin-flip contribution. Then we can see that the effect of 
the spin-flip contribution overwhelms the one of the spin-nonflip contribu- 
tion. The interaction contribution Xex always negative , and dominant over 
Xkin ^'^^ densities, while the kinetic contribution Xkm is always positive. 
If etc > 37r/2 = 4.7, X becomes negative over all densities. 
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For a given set of and Oc, x changes its sign at a certain density, 
denoted by rid, and it is a signal for the second-order phase transition. Note 
that the ferromagnetic transition in our case is of the first order, so that 
it is not sufficient to only see the magnetic susceptibility; even above that 
density the ferromagnetic phase may be possible. Actually there is a range, 
^ci < < nc2, where x > but e < 0. 

Above the density there is no longer the stable ferromagnetic phase. 
However, the metastable state is still possible up to the density rics, which 
is specified by the condition s.t. r] = detot/dp \p=i < 0. In Fig.2 we depict 
the quantities x~^:^^ and r] as the functions of density, e.g. for the set 
,mg = 300MeV and etc = 2.2. The crossing points with the horizontal axis 
indicate the critical deisities nci,nc2 and n^s, respectively. We can see that 
the ferromagnetic instability occurs at low densities, while the metastable 
state can exist up to rather high densities. 

Finally we show the critical lines satisfying = 0,6e = and r/ = in 
the QCD parameter and rriq) plane, which seperate the three character- 
istic regions for a given density. In Fig. 3 we demonstrate them at a density 
Uq = 0.3fm~^. All the lines have the maxima around the medium quark 
mass, and the mechanism of ferromagnetism is different for each side of the 
maximum, as already discussed. If we take niq = 300MeV for the s quark 
or TUq ~ OMeV for the u or d quark, and ac = 2.2 as in the MIT bag model 
again [22] , the quark liquid can be ferromagnetic as a metastable state. 

We have seen that the ferromagnetic phase is realized at low densities and 
the metastable state is plausible up to rather high densities for a reasonable 
range of the QCD parameters. Our calculation is based on the lowest-order 
perturbation. So we need to examine the higher-order gluon-exchange con- 
tributions to confirm the possibility. It should be interesting to refer a recent 
paper [SI], where the author also found the ferromagnetic transition at low 
densities within the perturbative QCD calculation beyond the lowest-order 
diagram. 

If a ferromagnetic quark liquid exists stably or metastably around or 
above nuclear density, it has some implications on the properties of strange 
quark stars and strange quark nuggets [SH]- They should be magnetized 
in a macroscopic scale. Considering a possibility to attribute magnetars 
to strange quark stars in a ferromagnetic phase, we roughly estimate the 
strength of the magnetic field at the surface of a strange quark star. Tak- 
ing the stellar parameters of strange quark stars to be similar to those for 
canonical neutron stars with the typical mass around Mq = IAMq, we find 
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Figure 3: Phase diagram in the couphng-strength (Q;s)-mass (mq) plane. 



the total magnetic dipole moment Mg, Mq = fig ■ (47r/3 • r^)nq for the quark 
sphere with the density Ug and the radius r^, where /ig is the magnetic mo- 
ment of each quark. Then the dipolar magnetic field at the star surface 
r = ~ 10km takes a maximal strength at the poles, 



3 \r) 



HgUg 
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15 



^O.lfm"^ 

with nuclear magneton which looks enough for magnetars. 



(16) 



2.2 Self-consistent calculation 

If we understand FM or magnetic properties of quark matter more deeply, 
we must proceeds to a self-consistent approach, like the Hartree-Fock theory, 
beyond the previous perturbative argument ^ . 
We begin with an OGE action: 



lint 2 / d ^ 



,(17) 



^Simple plane wave is the solution of the Hartree-Fock equation in the nonrelativistic 
electron gas, while it is not in quark matter. 
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where D'^'^ denotes the gluon propagator. By way of the mean-field approx- 
imation, we have 



hiF = J j^,Hp)G-/{pmp). (18) 

The inverse quark Green function G^^{p) involves various self-energy (mean- 
field) terms, of which we only keep the color singlet particle-hole mean-field 

V{p), 

Ga{p)-^ =i)-m + ii + V{p). (19) 

Taking into account the lowest diagram, we can then write down the self- 
consistent equations for the mean-field, V: 

-V{k) = {-rgf j ^^{-^D^'^{k-p)]^,^{-^GA{p)]^^^. (20) 

(A) 

Applying the Fierz transformation for the OGE action ()17|1 we can see 
that there appear the color-singlet scalar, pseudo-scalar, vector and axial- 
vector self-energies by the Fock exchange interaction. Taking the Feynman 
gauge for the gluon propagator, a manupilation gives 

{A) = ^^l-^^TTiGA)+^^,^T{GA^l,)-l[Y^TiGAl,)+^^^^ 

+ {color non— singlet or flavor non— singlet terms}. (21) 

When we restrict the ground state to be an eigenstate with respect to color 
and flavor, there is only left the flrst term which is color singlet and flavor 
singlet. Still we must take into account various mean-flelds in V , V = Us + 
ijbUps + 'jfj.UI^ + l/ybUa^ with the mean-flelds Ui. Here we only retain Ua^,(= 
U^) for simplicity and suppose that others to be vanished; 

r(A;)=775-UA(k), (22) 

with the static axial- vector mean-fleld ^/^(k). 

The poles of Ga{p), detG'^^(po + /^=en)=0, give the single-particle energy 
spectrum: 

en = ±e± (23) 



6±(p) = \V' + Ui(p) + m2 ± 2. /m2Ui(p) + (p ■ \J a{v))\ (24) 
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where the subscript in es(p), s = ± represents spin degrees of freedom, and 
the dissolution of the degeneracy corresponds to the exchange splitting of 
different "spin" states; the spectrum is reduced to a famihar form e± ~ 
m + ^ ± \ Va\ in the non-relativistic hmit [22] • 

There appear two Fermi seas with different volumes for a given quark 
number due to the exchange splitting in the energy spectrum. The appear- 
ance of the rotation symmetry breaking term, oc pUa in the energy spectrum 
implies deformation of the Fermi sea: thus rotation symmetry is violated in 
the momentum space as well as the coordinate space, 0(3) 0(2). Ac- 
cordingly the Fermi sea of majority quarks exhibits a "prolate" shape {F~), 
while that of minority quarks an "oblate" shape {F~^) as seen in Fig. 0] 






Pt Pz P: 



Figure 4: Modification of the Fermi sea as f/A(=const.) is increased from left 
to right. The larger Fermi sea {F~) takes a prolate shape, while the smaller 
one {F~^) an oblate shape for a given Ua- In the large Ua limit (completely 
polarized case), disappears as in the right panel. 

Then the self-consistent equation (j2(Jj) is reduced to the form. 



4iVc J (27r)^ ^ e,(p)^ - |k - e^(p) 

(25) 

with Pp = \/pl + rn? by taking along the z axis. Here we have discarded 
the contribution of the Dirac seas and only taken into account that of the 
Fermi seas. 

In the following we demonstrate some numerical results by replacing the 
original OGE by the "contact" (zero-range) interaction, D^^ —>■ —g'^'^/A'^, 
which may correspond to the Stoner model in the condensed matter physics 
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ESI- ^ 

We can easily see that the mean-field Ua becomes then momentum- 
independent, and the expression for Ua, Eq- fl^ - is proportional to the 
simple sum of the expectation value of the spin operator over the Fermi seas; 



c 



^tr7573G'A(p) 

d^p UA{p)+f3p 
p+ (27r)3 e+(p) 



+ 



(2vr)3 e_(p) 



.(26) 



2.3 Phase diagram on the temperature-density plane 

We will present the phase diagram in the three-fiavor case under two con- 
ditions jSZI: the chemical equilibrium condition (CEC) fiu = fJ'd = fJ^s and 
the charge neutral condition without electrons (CNC) Pu = Pd = ps, where 
quark masses are taken as rriu = Tn^ = 5MeV and = 150 — 350MeV, 

i.e., ps = -^J f^u d + "^s ~ ""^M d T = 0. In both conditions, since the spin 
polarization caused by the axial-vector mean-field is fully enhanced by the 
quark mass for given density or temperature, choice of the current quark 
mass seriously affects the results; especially, largeness of the strange quark 
mass has an essential effect on spin polarization. To get the phase diagram 
or critical line on the temperature-density plane, we use the thermodynamic 
potential Q within the mean-field approximation. 



n 



B=±l s=± i=u,d. 

d^k 



d^k 

(M 



-Tlog < exp 



es(k,mi, Ua) - Bpi 
T 



i=u,d,s 



(27r)3 



f/2 

e,(k,mi,[/^) + 



(27) 



where we have used the "contact" interaction, g = g /A , in place of the 
OGE interaction. Note that we take into accout the vacuum contribution 
in this formula (the second term in Eq. (jTTj) ). which should be regularized 
by ,e.g., the proper-time method (see §4.3). We can confirm that the ther- 
modynamic potential reproduces the self-consistent equation for the order 

^When we take into account the Debye screening, the time component of the gluon 
propagator becomes finite range due to the Debye screening. If typical momentum transfer 
Q is much smaller than the screening mass Ml, - Nfg'^^x^ I{2tt'^) gSI, we may replace the 
OGE with the infinite range by the zero-range effective interaction. 
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parameter Ua Eq. (j25j) in the three-flavor case, except the vacuum contri- 
bution. The vacuum (the Dirac sea) contribution always works against spin 
polarization as it should do, while the contribution of the Fermi sea gives 
rise to spontaneous spin polarization. 
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Figure 5: Phase diagram for two cases of strange quark masses, 250, 350MeV. 

Fig. El shows the critical temperature (the Curie temperature) as a func- 
tion of baryon-number density under the two conditions mentioned above; 
CEC and CNC. We can see that CNC tends to facilitate the system having 
spin polarization than CEC. This is because CNC holds the larger strange- 
quark density than CEC. 

Since the axial-vector mean-field arises from the Fock exchange inter- 
action among quarks in the Fermi sea and causes a kind of particle-hole 
condensation, there exists a critical density for a given coupling constant g. 
We show the critical density by varying the effective coupling constant g in 
Fig. ini The critical density is more lowered with the larger coupling strength, 
and this tendency is remarkable in the case of CNC. The result also indicates 
that even for the weak-coupling regime in QCD, spin polarization may appear 
at sufficiently large densities and low temperatures. 
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Comparison between CEC and CNC 




Figure 6: Critical density as a function of g under two conditions: CEC and 
CNC. 

3 Color magnetic superconductivity 
3.1 General framework 

If FM is realized in quark matter, it might be in the CSC phase. In this 
section we discuss a possibility of the coexistence of FM and CSC, which we 
call Color magnetic superconductivity [SB]. 

Recall the OGE action Eq. (fTTj) . By way of the mean-field approximation, 
we have 



IMF 



(27r)4 

in the Nambu-Gorkov formalism, allowing not only the particle-hole but also 
the particle-particle mean- field. The inverse quark Green function G~^{p) 
involves various self-energy (mean-field) terms, of which we only keep the 
color singlet particle-hole V{p) and color 3 particle-particle (A) mean-fields; 
the former is responsible to ferromagnetism, while the latter to superconduc- 
tivity. 



G-\p) 



:^-m + ^ + V{p) 7o A1" (p)7o_ 

A{p) f — m — ill + V{p) 

Gu{p) Gu{p) 

G2l{p) G22{P) 



-1 



(29) 
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where 

^^(k) = Cij^i-k), V = CV^C-\ (30) 

Taking into account the lowest diagram, we can then write down the self- 
consistent equations for the mean-fields, V and A: 

-V{k) = i-zgrj -^^{-^D^'^{k-p)h,^{-^Gn{p)h.Y■ (^l) 

and 

-m = {-W? j ^,{-^D^''{k-p)]^,^^{-^G,,{p)]^^^, (32) 
(c.f. Eq. (1201)). 

A = > /> % < ' < ^ 

G21 




Figure 7: Graphical interpretations of the coupled equations flHTj) and (jH^ 
with coefficients in front of R.H.S. given by N^.. The lower diagram becomes 
dominant in the large N^. limit. 

The structure of Eq. is the same as Eq. (pUj) . and we can see there 
appear the color-singlet scalar, pseudoscalar, vector and axial-vector self- 
energies by applying the Fierz transformation . Here we retain only Us,U^, U^^ 
in V and suppose that others to be vanished as before. We shall see this 
ansatz gives self-consistent solutions for Eq. ()31|) within the zero-range ap- 
proximation for the OGE interaction because of axial and reflection symme- 
tries of the Fermi seas. We furthermore discard the scalar mean-field Ug and 
the time component of the vector mean-field f/° for simplicity since they are 
irrelevant for the spin degree of freedom. 

According to the above assumptions and considerations the mean-field V 
in Eq . fl29|) renders 

^ = 7375^7^, Ua^UI, (33) 
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with the axial- vector mean- field Ua, as in Eq. (j22|) . Then the diagonal com- 
ponent of the Green function G'ii(p) is written as 



Gnip) 



-1 



(34) 



with 



(35) 
(36) 



where 7573 = 7573 and Ga{p) is the Green function with Ua which is deter- 
mined self-consistently by way of Eq. (^0]) . 



3.2 type anisotropic pairing 

Before constructing the gap function A, we first find the single-particle spec- 
trum and their eigenspinors in the absence of A, which is achieved by diag- 
onalization of the operator G^. We have already known four single-particle 
energies t± (positive energies) and —e± (negative energies), which are given 
as 

e±(p) = + ^72 + ^2 ± 2UA^f^^^^T^,, (37) 

and the eigenspinors 0s, s = ± should satisfy the equation, (^^^(es, p)0s = 0. 
Here we take the following ansatz for A: 

A(p) = 5^A,(p)5,(p), 

s=± 

Bsip) = lo<P-s{p)<Pl{p). (38) 

The structure of the gap function (j38|) is then inspired by a physical 
consideration of a quark pair as in the usual BCS theory: we consider here 
the quark pair on each Fermi surface with opposite momenta, p and — p so 
that they result in a linear combination of J'" = 0~, 1~ (see Fig. |H1) ^ • 

^Recently spin-one color superconductivity has been also studied in the normal matter 

m 
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Figure 8: Deformed Fermi seas and 
the quark pair on each surface. The 
top figures show those in the absence 
of A± and the middle figures diffusion 
of the Fermi surfaces in the presence 
of A±. The bottom ones show the 
quark pairing on the Fermi surfaces. 



Figure 9: Schematic view of the polar- 
angle dependence of the gap functions 
at the Fermi surface, (a) for m = 
and (b) for m 7^ 0. 



As is still a matrix in the color-flavor space. Since the anti-symmetric 
nature of the fermion self-energy imposes a constraint on the gap function 

0, 

CAip)C-' = A^(-p). (39) 

An{p) must be a symmetric matrix in the spaces of internal degrees of free- 
dom. Taking into account the property that the most attractive channel of 
the OGE interaction is the color anti-symmetric 3 state, it must be in the 
flavor singlet state. Thus we can choose the form of the gap function as 

(A,) = e'^^h'^As (40) 
for the two-flavor case (2SC), where a, (3 denote the color indices and i, j the 
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flavor indices. Then the quasi-particle spectrum can be obtained by looking 
for poles of the diagonal Green function, Gn: 

p _ / v/ (6.(p)-/xP + ^:W for color 1, 2 

I for color 3 ^ ^ 

Note that the quasi-particle energy is independent of color and flavor in this 
case, since we have assumed a singlet pair in flavor and color. 

Gathering all these stuffs to put them in the self- consistent equations, we 
have the coupled gap equations for A^, 

AAk,e,) = ^9 J^p sm9pJ2Ts'sik,9,,p,ep)^^^^, (42) 
and the equation for Ua, 

within the "contact" interaction, = g'^/A'^, where ff(p) denotes the mo- 
mentum distribution of the quasi-particles. We find that the expression for 
Ua, Eq- (j^Hj) . is nothing but the simple sum of the expectation value of the 
spin operator with the weight of the occupation probability of the quasi- 
particles Vg for two colors and the step function for remaining one color (cf. 

Carefully analyzing the structure of the function Tg's in Eq. (j42j) . we 
can easily find that the gap function should have the polar angle [6) 
dependence on the Fermi surface. 



R + F 



^ym"^ + {pf{0) COS 



(44) 



with constants F and R to be determined (see Fig. 

As a characteristic feature, both the gap functions have nodes at poles 
{6 = 0, vr) and take the maximal values at the vicinity of equator {6 = 
7r/2), keeping the relation, A_ > A+. This feature is very similar to 
pairing in liquid ^He or nuclear matter jJOJ E]; actually we can see our 
pairing function Eq. (jl^ to exhibit an effective P wave nature by a genuine 
relativistic effect by the Dirac spinors Accordingly the quasi-particle 

distribution is diffused (see Fig. ^ 
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We demonstrate some self-consistent solutions here. Since we have little 
information to determine the values of the parameters g and 6 (there may 
be other more reasonable form factors than the present cut-off function), 
and our purpose is to figure out qualitative properties of spin polarization in 
the color superconducting phase, we mainly set in the following calculations 
them as ^ = 0.13 MeV~^ and 6 = 0.1/^, for example, which is not so far from 
the couplings in NJL-like models [HI HSl EH • 

We first examine spin polarization in the absence of CSC. In Fig. ITUl we 
show the the axial-vector mean-field Ua, with A± being set to be zero, as a 
function of baryon number density Pb{= Pql'i) relative to the normal nuclear 
density po = 0.16 fm~'^ for m = 14 ~ 25 MeV (dashed lines). It is seen 
that the axial-vector mean-field (spin polarization) appears above a critical 
density and becomes larger as baryon number density gets higher. Moreover, 
the results for different values of the quark mass show that spin polarization 
grows more for the larger quark mass. This is because a large quark mass 
gives rise to much difference in the Fermi seas of two opposite "spin" states, 
which leads to growth of the exchange energy in the axial- vector channel. 

Next we solve the coupled equations (ji^ and with Eq. (jH)). Results 
for Uai R and F are shown in Fig. (solid lines) 

As a consequence, we can say that FM and CSC barely interfere with 
each other 

3.3 Another possibility - Gapless type pairing 

Nowadays there have been many studies about the pairing of quarks in the 
two Fermi spheres with different sizes, which is caused by the mass and charge 
differences among three-flavor quarks. It is well known that fermion pairing 
between two different Fermi surfaces gives rise to the LOFF phase j45[ HB] or 
the gapless superconducting phase EZl HH] • There have been discussions 
about the phase separation and the mixed phase in these context 130]. 

In the presence of magnetization we have seen that there are two Fermi 
seas with different size and deformation, depending on the spin polarization. 
So we can consider another pairing than the previous one: two quarks with 
opposite momenta and polarizations with each other take part in the pairing 
jSlj . Introducing the following notation, 

e„ = {e„,e+,-e_e+} (n = 1 ~ 4) (45) 
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Figure 10: Axial-vector mean-field as a function of baryon number density 
(po = 0.16 fm-3) for g = 0.13 MeV'^ and 5 = 0.1/i. (a) for m = 14 ~ 16 
MeV and (b) for m = 20 and 25 MeV. Dashed (Solid) lines are obtained in 
the normal (color superconducting) phase. 

for the single quark energy by using Eq. (j37j) . The pairing function can be 
written in the similar form to Eq. (j38|) 

4 

A(p) = ^A„(p)5„(p) (46) 

n=l 

with 

i^n(p)=7o0-n(p)0|,(p), (47) 

where 0_fj(p) is defined by 0_^j(p) = 0-i+(-i)"(p)- As a combination of the 
quark pair in the color and flavor spaces, we assume it to be anti-symmetric 
in both spaces, 

[MP)U■^j = ea(33eijA„{p) (48) 

as before. Then the quasi-particle energy is given by 



p .„^ _> ±V(i^„^)^+ |An(p)P for color 1, 2 
^-^^^^ - ^ ±V (e„,(p)-/.)^ for color 3 ^^^^ 
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with 

ES,a ^ {e^-f^)T{e,-^,) ^ ^^^^ 

which clearly exhibits a gapless excitation. We can also see that the gap 
function shows cos6'-like dependence on the Cooper surface defined by the 
equation, E^{p, 0) = 0. These features resemble those given in ref. [S2], where 
the electron pairing with spin anti-parallel component of the S = 1 triplet is 
considered in the presence of magnetization. 



4 Dual chiral density wave 

4.1 Chiral symmetry restoration and Instability of the 
directional mode 

We consider here another type of magnetism in quark matter at moderate 
densities, which is closely connected with chiral symmetry. We shall see 
that the ground state in the spontaneously symmetry breaking (SSB) phase 
becomes unstable with respect to producing a density wave. Accordingly 
the quark magnetic moment spatially oscillates and a kind of spin density 
wave is induced. The density wave can be described as a dual standing wave 
in the scalar and pseudo-scalar densities where they spatially oscillate 
in the phase difference of 7r/2 to each other. It is well known that chiral 
symmetry is spontaneously broken due to the quark (g)-anti-quark (g) pair 
condensate in the vacuum and at low densities; since we take the vacuum as 
an eigenstate of parity operation, only the scalar density is non vanishing to 
generate finite mass of quarks. Geometrically both the scalar and pseudo- 
scalar densities always reside on the chiral sphere with the finite modulus in 
the SSB phase, and any chiral transformation with a constant chiral angle 
da shifts each value on the sphere, leaving the QCD Lagrangian invariant. 
The spatially variant chiral angle 6{r) represents the degree of freedom of the 
Nambu-Goldstone mode in the SSB vacuum. The dual chiral density wave 
(DCDW) is described by such a chiral angle 0{r). When the chiral angle has 
some space-time dependence, there should appear extra terms in the effective 
potential as a consequence of chiral symmetry: one trivial term is the one 
describing the quark and DCDW interaction due to the non-commutability 
of ^(r) with the kinetic (differential) operator in the Dirac operator. Another 
one is nontrivial and comes from the vacuum polarization effect: the energy 
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spectrum of the quark is modified in the presence of 6{r) and thereby the 
vacuum energy has an additional term, oc (V0)^ in the lowest order. This 
can be regarded as an appearance of the kinetic term for DCDW through the 
vacuum polarization [31]. Thus, the interaction becomes strong enough to 
overwhelm the kinetic energy increase, the state becomes unstable to generate 
DCDW. 

Many studies have suggested that chiral symmetry is restored at a certain 
density by suppression of qq excitation due to the presence of the Fermi 
sea, where none of the mean- fields is present. In usual discussion of such 
symmetry restoration, one implicitly discards the pseudo-scalar mean-field 
and is concentrated in the behavior of the scalar mean-field, while there is no 
compelling reason for the pseudo-scalar density to be vanished. Allowance of 
the degree of freedom of the chiral angle is nothing else but the appearance of 
DCDW. Thus we can say that instability of the ground state with respect to 
forming DCDW provides another path to symmetry restoration (see Fig.lll|). 




Figure 11: Schematic view of the possible paths of symmetry restoration in 
the chiral space. There may be possible another paths utilizing the degree 
of freedom of the chiral angle, besides the usual one along 6 = 0. 
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4.2 DCDW in the NJL model 



Taking the Nambu-Jona-Lasinio (NJL) model as a simple but nontrivial ex- 
ample, we explicitly demonstrate that quark matter becomes unstable for a 
formation of DCDW above a certain density; the NJL model has been re- 
cently used as an effective model of QCD, embodying spontaneous breaking 
of chiral symmetry in terms of quark degree of freedom j55j ^ . We shall 
explicitly see the DCDW state exhibits a ferromagnetic property. 

We start with the NJL Lagrangian with Nf = 2 flavors and Nc = 3 colors, 

Cnjl = - m,)i; + G[{^Pijf + (iji-f^r^jf], (51) 

where rric is the current mass, rric — 5MeV. Under the Hartree approximation, 
we linearize Eq. ()51|) by partially replacing the bilinear quark fields by their 
expectation values with respect to the ground state. 

In the usual treatment to study the restoration of chiral symmetry at finite 
density, authors implicitly discarded the pseudo-scalar mean-field, while this 
is justified only for the vacuum of a definite parity. We assume here the 
following mean-fields, 

(ipip) = Acos(q-r) 
(^^75T"3^) = Asin(q-r), (52) 

and others vanish . This configuration looks to break the translational in- 
variance as well as rotation symmetry, but the former invariance is recovered 
by absorbing an additional constant by a global chiral transformation. Ac- 
cordingly, we define a new quark field ipw by the Weinberg transformation 

ipw = exp[i75r3q ■ r/2]?/', (53) 

to separate the degrees of freedom of the amplitude and phase of DCDW in 
the Lagrangian. In terms of the new field the effective Lagrangian renders 

Cmf = ^PwW -M- l/2-f,T3i]^w - GA\ (54) 

^We can see that the OGE interaction gives the same form after the Fierz transforma- 
tion in the zero-range limit 38] 

^It would be interesting to see that the DCDW configuration is similar to pion condensa- 
tion in high-density nuclear matter within the a model, considered by Dautry and Nyman 
fDN) |5(jl 1571 158| . where a and tt^ meson condensates take the same form as Eq. (j52|l . The 
same configuration has been also assumed for non-uniform chiral phase in hadron matter by 
the use of the Nambu-Jona-Lasinio model [SUj. However, DCDW is by no means the pion 
condensation but should be directly considered as particle-hole and particle-antiparticle 
quark condensation in the deconfinement phase. 
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where we put M = —2GA and q'^ = (0, q), taking the chiral hmit (m^ = 0). 
The form given in (j54p appears to be the same as the usual one, except the 
axial-vector field generated by the wave vector of DCDW; the amplitude of 
DCDW produces the dynamical quark mass in this case. We shall see the 
wave vector q is related to the magnetization: the phase of DCDW induces 
the magnetization. With this form we can find a spatially uniform solution 
for the quark wave function (see Table Ej), ipw = uw{p) exp(zp ■ r), with the 
eigenvalues, 

= ^JeI + |q|74 ± v/(p ■ q)2 + E, = (M^ + (55) 

for positive-energy (valence) quarks with different spin polarizations (c.f. 

dSII)). ® 

Table 3: Diagram of the Weinberg transformation. 

(^^) ^ ^ {^w^w) = A(^ 0) 

{ijji'j5'T'3i') ^ {ipwi^bTsipw) = 

q/2ocVe ("axial- vector") 
non-uniform uniform 



4.3 Thermodynamic potential 

The thermodynamic potential is given as 

fitotai = 7 / ^ E [(^P - /^)^^ - + ^V4G 

= + ^^vac + MV4G. (56) 

where e± = 6(12 - E^), fx is the chemical potential and 7 the degeneracy 
factor 7 = NfNc- The first term f2vai is the contribution by the valence 
quarks filled up to the chemical potential, while the second term f2vac is 
the vacuum contribution that is apparently divergent. We shall see both 

®This feature is very different from refs.|16|. where wave function is no more uniform. 
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contributions are indispensable in our discussion. Once fi^otai is properly 
evaluated, the equations to be solved to determine the optimal values of A 
and q are 

^ = ^ = °- 

Since NJL model is not renormalizable, we need some regular izat ion pro- 
cedure to get a meaningful finite value for the vacuum contribution Q^ac, 
which can be recast in the form, 

^^vac = y J2^^^^^^^' ^^^^ 

with use of the propagator Sw = {][> — M — 1/2x375;^)^^. There are various 
kinds of regularization and we must carefully choose the relevant one to the 
theoretical framework. Since the energy spectrum is no more rotation sym- 
metric, we cannot apply the usual energy or momentum cut-off regularization 
(MCOR) scheme to regularize VL^^c- Moreover, the regularization should be, 
at least, independent of the order parameters A and q. Note that this de- 
mand is essential to discuss the phase transition: improper regularizations 
spoil the consistency of the framework and give unphysical results for the 
order parameters A and q through Eq. (jSTjl . We adopt here the proper-time 
regularization (PTR) scheme [HII, which is one of regularizations compatible 
with Eq. ()57|1 ^ . Introducing the proper-time variable r, we eventually find 



J r dr r dpz 



^vac 



87^3/2 / ^5/2 /_ 271 



~(^/^|TM^+g/2)2r ^ g-(^/^|^-g/2) Vj 

(59) 

which is reduced to the standard formula in the limit g — 0. 

The integral with respect to the proper time r is still divergent due to the 
r ~ contribution. Regularization proceeds by replacing the lower bound of 
the integration range by 1/A^, which corresponds to the momentum cut-off 
in the MCOR scheme. 

Now we examine a possible instability of quark matter with respect to 
formation of DCDW. In the following we first inquire the sign change of the 
curvature of fitotai at the origin {stiffness parameter), [3. Expanding f^^ac 
with respect to q up to O(g^), we find 

n,ac = ^^°ac + /5vacg' + 0(g^) (60) 



^The Pauli-Villars reguralization may be another candidate. 
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where the vacuum stiffness parameter /^vac is given by 



/3vac = -^J(MVA') (61) 

with a universal function, J{x) = — xEi(— x). The nontrivial term originates 
from a vacuum polarization effect in the presence of DCDW and provides 
a kinetic term (oc (VOy) for DCDW . The vacuum stiffness parameter /5vac 
can be also written as = \fl jSHj with the pion decay constant and 
is always positive; it gives a 'repulsive' contribution, so that the vacuum is 
stable against formation of DCDW. Note that it gives a null contribution in 
case of M = , irrespective of g, as it should be. 

For given /i, M and q we can evaluate the contribution by the Fermi seas 
i^vai using Eq. but its general formula is very complicated jH31- However, 
it may be sufficient to consider the small q case for our present purpose. Then 
the thermodynamic potential can be expressed as 

^vai = nl,,-^,M\^H{y^/M) + 0{q') 

= ^L + firar + 0(g^) (62) 

up to 0{q^), where H{x) = \ii{x + \/x^ — 1) and = Cvai^/^Pvai "^i^^ p°^i = 
-^{fJ''^ — IVPY^"^ for normal quark matter. The valence stiffness parameter 
then reads 

f^vai = -^M^H{i^/M) (63) 

Since the function H{x) is always positive and accordingly P^ai < 0, the 
magnetic term fi™^^ always gives a negative energy and approaches to zero 
as M ^ (triviality). 

We may easily understand why the valence quarks always favor the for- 
mation of DCDW. First, consider the energy spectra for massless quarks (see 
Fig.d. 

As is already discussed, our theory becomes trivial in this case and we 
find two spectra 

E^ = \/pl + i\Pz\±<im', P± = iPx,Py,0), (64) 

which are essentially equivalent to = |p| with definite chirality. 

There is a level crossing at p = 0. Once the mass term is taken into 
account this degeneracy is resolved and the energy splitting arises there. 
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Figure 12: Energy spectra for p_|_ = 0. with M = (thick sohd and 
dashed hnes). with the definite chirahty is also shown for comparison 
(dotted fine). We can see there is a degeneracy of E"^ at Pz = for M = 0, 
while it is resolved by the mass (thin solid and dashed lines). 

Hence it causes an energy gain, if g = 0(2yu); we can see that this mechanism 
is very similar to that of SDW by Overhauser 120] . 

Using Eqs. (fHUI). (jH^ we write the thermodynamic potential as 

^^total = ^^7VJL + /5g' + 0(g^) (65) 

with the total stiffness parameter j3 = Pvac+Pvai and the usual NJL expression 
without DCDW, n^jL = ^^°^,(M) + ^]0^l(M) + MV4G. The dynamical quark 
mass M is given by the equation, Sfitotai/f?^ = 0; At the order of g° the 
dynamical quark mass M° is determined by the equation, dflNji/dMlj^p = 
0. Since M — M° = 0{q'^), DCDW onsets at a certain density where the total 
stiffness parameter f3 becomes negative: the critical chemical potential fi'^^ is 
determined by the equation, 

(3 = \fl - ^ {M'f H{f^-/M') = 0. (66) 

Note that this is only a sufficient condition for formation of DCDW, and 
we can never exclude the possibility of the first-order phase transition or 
metamagnetism [TTJ 122. Actually, we shall see that DCDW occurs as a 
first-order phase transition. 
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Figure 13: Contours of fitotai at T = are shown in M — g plane as the 
chemical potential increases, (a) — > (f). The cross in each figure denotes the 
absolute minimum. 



4.4 First-order phase transition 

The values of the order parameters M and q are obtained from the minimum 
of the thermodynamic potential (j56j) for T = 0. Fig. IT^ shows the contours 
of fitotai in the M-q plane as the chemical potential increases, where the 
parameters are chosen as GA^ = 6 and A = 850 MeV, to reproduce the 
constituent quark mass in the vacuum {ji = 0) [53] . 

The crossed points denote the absolute minima. There are two critical 
chemical potential /x = fici, f^c2- for the lower densities (Fig. [Tnfa)-(b)) the 
absolute minimum resides at the point (M 7^ 0, g = 0) indicating the SSB 
phase. At fi = /id (Fig. [THfc)) the potential has the two absolute minima at 
(M 7^ 0, g = 0) and (M 7^ 0, g 7^ 0), showing the first-order transition to the 
DCDW phase which is stable for fid < fi < fic2 (Fig. (fT^ d-e). At /i = /ic2 
(Fig. IT^ f)) the axis of M = and a point (M 7^ 0, g 7^ 0) become minima, 
the system undergoes the first-order transition again to the chiral-symmetric 
phase. 

Fig. summarizes the behaviors of the order-parameters M and g as 
functions of /i at T = 0, where that of M without DCDW is also shown for 
comparison. It is found that DCDW develops at finite range of fi (/id < 
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^ f^c2), where the wave number q increases with /i but its value is smaller 
than twice of the Fermi momentum 2/c^(~ 2^ for free quarks) since the 
nesting of Fermi surfaces is incomplete in the present 3-D system; actually, 
the ratio becomes q/ kp = 1.17—1.47 for the baryon-number densities pb/po = 
3.62 - 5.30 where DCDW is stable (see Fig.CH). 

Order Parameters at 7=0 
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Figure 14: Wave number q and the 
dynamical mass M are plotted as 
functions of the chemical potential 
at T = 0. Solid (dotted) line for M 
with (without) the density wave, and 
dashed line for q. 
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Figure 15: Baryon number density as 
a function of p. po = 0.16fm~^: the 
normal nuclear density. 



4.5 Correlation functions 

In this section, we consider scalar- and pseudoscalar-correlation functions, 
ns,sp(^), in the massless limit M — > 0, and discuss their relation with the 
mechanism for DCDW. In the static limit ko —>■ 0, the correlation functions 
have a physical correspondence to the static susceptibility for the spin- or 
charge- density wave ^Tj. We shall see that these functions have a differential 
singularity at k = 2k p, reflecting the sharp Fermi surface at T = 0. 

We explicitly evaluate the effective interactions, Tg^gp^k), in the pseudo- 
scalar and scalar channels within the random phase approximation ^1 , 
which are related to the correlation functions ns,sp(fc), i.e., 2GIls,sp{k) = 
^ s,sp{k)Ilssp{k): 

^rs,ps(fc) = ,_,l%jky ^'^^ 
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where ngpg(fc) are the polarization functions in medium, 



2kl)-2NfN,^^I{k^)\M^o 

in the static and chiral hmit [^2] • It is well known that poles of the effective 
interaction give the energies of scalar and pseudo-scalar mesons jlH ESI • Note 
that the inverse of the effective interaction in the massless limit also gives 
the coefficient of in the effective potential in the presence of DCDW, 

^otal = ^otallM^O + ^T^,' iq)\M^oM^ + 0{M'). (69) 

Hence, the critical density and the critical wave vector gcrit should be given 
by the equations, 

r-i(fet)|Af^o = 0, 5rp,i(g)|M^o/5g|... = o, (70) 

in the case of the second-order or weakly first-order phase transition. FiglT^ 
show the function l/rps(|k|)|jvf=o and we can see these conditions are almost 
satisfied at the terminal density, /x = /ic2, due to a tiny jump in the dynamical 
mass (weakly first-order phase transition) . A numerical calculation in the 
chiral limit gives fit{= kp) = 0.5320A and |kt| = lAdSkp, which almost 
coincide with the previous results given in Figs. IT^ and fT^ yUc2 = 0.53254A 
and q = 1.469^^7 (where kp = M = 0.034A). On the other 

hand, the phase transition is of first order at the onset density, /i = /id: there 
is a discontinuous jump in the dynamical mass (the amplitude of DCDW) 
M, so that the above argument cannot be applied any more. Besides, the 
correlation functions or the effective interactions provide a powerful tool to 
analyze the DCDW phase as far as the dynamically generated mass (the 
amplitude of DCDW) M is small, as in the present case. From the behavior 
of the function rps(|k|)~^ shown in Fig. UHl it is found that rps(|k|)"^ takes 
the lowest value at |k| ~ 1.3 — 1.5kp{0{2kp)), reflecting the sharp Fermi 
surface, and thus a finite wave number q gives the lower potential energy 

""^^ Actually the function rps(|k|)~^ diverges at |k| = 2kp in the one-dimensional case, 
which means the complete nesting. 
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(a) m^=0 



(b) m<-=10 MeV 




Figure 16: Function, l/rps(|k|), is plotted for various Fermi momenta, 
kp/A = 0.4752,0.5, and 0.6. The thick (thin) hues correspond to the chiral 
hmit rric = (m^ = 5MeV). In the case of kp/A > 0.4752, the mass-gap 
equation has a extremum solution M = in the absence of DCDW. 

in Eq. (jU^ than g = in the density range of DCDW. These values are 
consistent with those in Fig. El and we can see again that DCDW is closely 
related to the sharpness of the Fermi surface. 

It should be noted again that the negative value of the function rps(|k|)~^| 
gives a necessary condition for formation of DCDW, but the sign change does 
not necessarily imply the critical condition in the case of first-order phase 
transitions, as in the present case; the terminal transition is weakly first- 
order and we can apply there. It should be also noted that its minimum 
point always gives an optimal value of the wave vector in the presence of 
DCDW. Thus we can see by the use of the correlation functions that the 
particle-hole pairing with finite momentum q = 0{2kF) effectively lowers the 
free energy in comparison with the zero total momentum. 

The above argument might also be available even for the case of a fi- 
nite current-quark mass, ~ 5MeV: Fig. shows that the minimum of 
Fpsdkl)"^ has little shift from that in the chiral limit. 



34 



4.6 Magnetic properties 

The mean-value of the spin operator is given by 

^ 1 t ^ lg/2±^ 

= -ul^E^uw = 2 — ^± ^ vac, (71) 

with l3.p = y/pf+lTT?, where "vac" means the vacuum contribution. First 
note that the integral of over the Fermi seas should be proportional to 
q, and the solution with q ^ seems to imply FM. However, we can show 
that PTR gives the vacuum (the Dirac sea) contribution oppositely to cancel 
the total mean- value of the spin operator, which is consistent with Eq. ()57|) . 
Instead we can see that the magnetization spatially oscillates, 

= {qauq) = {jqctu) cos(q ■ r), (72) 

with 

Jf+-f- {^T^r +pI 
which means a kind of spin density wave jST| . 

4.7 Phase diagram in the T — fi plane 

To establish the phase diagram in the T — /i plane, we derive the thermo- 
dynamic potential at finite temperature in the Matsubara formalism. The 
partition function for the mean-field Hamiltonian is given by 



J DipDijjexp J drj d^r i^tp i(9 + Mexp (i75q ■ r)— 70/i 
= n {(-n + ^)^-i^.^(k)}^^^^xexp|-(^)l//3|, (74) 

where P = 1/T, d = — 7o<9r + 27 V and Un the Matsubara frequency. Thus 
the thermodynamic potential is obtained, 

n0{q,M) = ~T log Zp{q,M)/V 

= -iV^iV,y"^5^{T log [e-^(i^^(k)-M)+i] [e-^(^^«+^)+l]+E,(k)} 
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From the absolute minima of the thermodynamic potential (j75|) . it is 
found that the order parameters at T 7^ behave similarly to those at T = 
as a function of n, while the chemical-potential range of the DCDW at finite 
temperature, fJ^ciiT) < /i < fic2iT), gets smaller as T increases. We show the 
resultant phase diagram in Fig. El where the ordinary chiral-transition line 
is also given. Comparing phase diagrams with and without q, we find that 
the DCDW phase emerges in the area (closed area in Fig. I17|) which lies just 
outside the boundary of the ordinary chiral transition. We thus conclude 
that the DCDW is induced by finite-density contributions, and has an effect 
to expand the chiral-condensed phase (M 7^ 0) toward low temperature and 
high density region. 



0.20 
0.15 

< 

0.10 
0.05 
0.00 

0.40 0.45 0.50 

Figure 17: Phase diagram obtained from the thermodynamic potential (j75|) . 
The solid (dashed) line shows the chiral-restoration boundary in the presence 
(absence) of DCDW. The closed area denoted by "DCDW" shows the DCDW 
phase. 



5 Summary and Concluding remarks 

We have seen some magnetic aspects of quark matter: ferromagnetism at high 
densities and spin density wave at moderate densities within the zero-range 
approximation for the interaction vertex. These look to follow the similar 
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development about itinerant electrons: Blocli mechanism at low densities 
and spin density wave at high densities by Overhauser. 

By a perturbative calculation with the OGE interaction, we have seen 
ferromagnetism in quark matter at low densities (§2.1) . It would be worth 
mentioning that another study with higher-order diagrams qualitatively sup- 
ports it 34 . These studies suggest an opposite tendency to the one using 
the zero-range interaction (§§2.2,2.3). Note that we can also see the same 
situation for itinerant electrons; the Hartree-Fock calculation based on the 
infinite-range Coulomb interaction favors ferromagnetism at a low density 
region, while the Stoner model, which introduces the zero-range effective in- 
teraction instead of the Coulomb interaction, gives ferromagnetism at high 
densities. So we must carefully examine the possibility of ferromagnetism in 
quark matter by taking into account the finite-range effect. 

We have examined the coexistence of spin polarization and color super- 
conductivity by choosing a quark pair with the same polarization. We have 
introduced the axial-vector self-energy and the quark pair field (the gap func- 
tion), whose forms are derived from the one-gluon-exchange interaction by 
way of the Fierz transformation under the zero-range approximation. Within 
the relativistic Hartree-Fock framework we have evaluated their magnitudes 
in a self-consistent manner by way of the coupled Schwinger-Dyson equations. 

As a result of numerical calculations spontaneous spin polarization occurs 
at a high density for a finite quark mass in the absence of CSC, while it never 
appears for massless quarks as an analytical result. In the spin-polarized 
phase the single-particle energies corresponding to spin degrees of freedom, 
which are degenerate in the non-interacting system, are split by the exchange 
energy in the axial-vector channel. Each Fermi sea of the single-particle 
energy deforms in a different way, which causes an asymmetry in the two 
Fermi seas and then induces the axial- vector mean- field in a self-consistent 
manner. In the superconducting phase, however, spin polarization is slightly 
reduced by the pairing effect; it is caused by competition between reduction 
of the deformation and enhancement of the difference in the phase spaces of 
opposite "spin" states due to the anisotropic diffuseness in the momentum 
distribution. 

We have also noted another possibility of the pairing: the quark pair with 
opposite polarization to each other. It may lead to a gapless superconductor, 
but we need a further study. 

We have seen that dual chiral desnity wave (DCDW) appears at a certain 
density and develops at moderate densities (§3). It occurs as a result of the 
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interplay between the qq and particle-hole correlations. The phase transition 
is of weakly first order, and the restoration of chiral symmetry is delayed 
compared with the usual scenario. 

For the discussion of DCDW given in §4, we have seen the remarkable 
roles of the Fermi sea and the Dirac sea: the former always favors DCDW, 
while the latter works against it. The similar situation also appears about the 
magnetic property of quark matter. The mean value of the spin operator over 
the Fermi seas of valance quarks always gives a finite value in the presence 
of DCDW, which is a kind of ferromagnetism, but the vacuum contribution 
given by the Dirac seas completely cancels it. As a result there is no net spin 
polarization in this case, but we have seen magnetization spatially oscillates 
instead (spin density wave) . This is one of the typical examples in which the 
nonrelativistic picture is qualitatively different from the relativistic one by 
the vacuum effect. 

It would be interesting to recall that DCDW is similar to pion conden- 
sation within the a model, considered by Dautry and Nyman jJHl! where cr 
and tt" meson condensates take the same form as Eq. (j52p . So it might be 
intriguing to connect pion condensation before deconfinement with DCDW 
after it in light of symmetry consideration. Note that this type of hadron- 
quark continuity has been also suggested in the context of hadron and quark 
superconductivities [S^ . 

If ferromagnetism is realized in quark matter, it may give a microscopic 
origin of the magnetic field in compact stars; actually we have seen that it 
can give a possible explanation for the superstrong magnetic field observed 
in magnetars, if they are quarks stars. It would be challenging to explain 
other characteristic phenomena in magnetars such as a sudden braking down 
observed in a soft gamma-ray repeater SGR 1806-20 or SGR 1900+14 
some global reconfiguration of of the magnetic has been suggested for these 
phenomena |64|. It would be also ambitious to give a scenario based on 
magnetic properties of quark matter, which can explain the hierarchy of 
the magnetic field observed in three classes of neutron stars, magnetars, 
radio pulsars and recycled millisecond pulsars. Ferromagnetism may give a 
permanent magnetization and there is no field decay, in difference from the 
dynamo mechanism caused by the charged current. 

The magnetic phases considered here accompany the symmetry breaking, 
5*0(3) 0(2), so that we can expect the Nambu-Goldstone modes as lowest 
excitations in the ground state: spin wave in ferromagnetism and phason in 
DCDW. It should be interesting to study these modes. Such low excitation 
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modes may affect the thermal evolution of compact stars jHSl- It would be 
also interesting to investigate how the effective interaction by exchanging 
such excitations between quarks affects superconductivity. 
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